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Abstract 

Stellar multiplicity is an ubiquitous outcome of the star formation process. Characteriz- 
ing the frequency and main characteristics of multiple systems and their dependencies on 
primary mass and environment is therefore a powerful tool to probe this process. While 
early attempts were fraught with selection biases and limited completeness, instrumenta- 
tion breakthroughs in the last two decades now enable robust analyses. In this review, we 
summarize our current empirical knowledge of stellar multiplicity for Main Sequence stars 
and brown dwarfs, as well as among populations of Pre-Main Sequence stars and embedded 
protostars. Clear trends as a function of both primary mass and stellar evolutionary stage 
are identified that will serve as a comparison basis for numerical and analytical models of 
star formation. 

1 INTRODUCTION 

The existence of binary and multip le systems was recognized soon after the birth of 
modern astronomy in the 17 th century. iKuiperl ( 1935aT ) first argued that empirically deter- 



mining the multiplicity frequency and distribution of key orbital parameters would prove 
to be highly valuable from a theoretical standpoint. The dependency of the multiplicity 
frequency and their associated orbital parameters on primary mass should contain the im- 
print of the physical processes at play throughout the lifetime of stellar populations. In 
particular, if dynamical interactions (within multiple systems, in stellar clusters and in the 
Galactic tidal field) can be properly accounted for, statistical properties of multiple sys- 
tems can serve as a posteriori tests for star formation theories. For instance, whether star 
formation is characterized by universal multiplicity distributions (akin to the apparently- 
universal Initial Mass Function, IMF, in the present day Universe) is an open question that 
has important bearing on star formation theory. 

The frequency and properties of multiple systems have been discussed in mul tiple entries 



in t he Annual Review : van de KamrJ ( 1975 ) discussed astrometric binaries, while lAbtl (| 19831 ) 



and iMathieul (|l994T ) presented statistical analyses of multiplicity among Main Sequence 



(MS) field stars in and Pre-Main Sequence (PMS) stars, respectively. Knowledge of stellar 
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multiplicity has greatly expanded since these seminal reviews, due to improvements in 
instrumentation and larger, more uniform s amples. In addition, new stellar populations 
have been probed, such as high-mass stars Jzinnecker fc Yorkel 2007 ) and very low-mass 
(VLM) stars and brown dwarfs (BDs, see Luhmanl 2012 ). allowing for a more complete 
picture to be drawn. On the theoretical side, prompt fragmentation during the gravitational 
collapse of a prestellar core appears to be the favored formation mecha nism although no 
single process appears to account for all observed multiplicity properties (|Tohlinell2002H . 

This review summarizes our current empirical knowledge of the frequency and principal 
properties (such as the distributions of orbital periods, mass ratios, and eccentricities) of 
multiple systems as a function of stellar mass and age. Multiple systems are thought to 
play central roles in numerous aspects of post-MS evolution, such as the blue straggler 
phenomenon, X-ray binaries, Type la supernovae, and the shaping of bipolar planetary 
nebulae. However, we restrict the scope of this review to MS and less evolved objects, 
where unbiased uniform surveys have targeted larger samples. 



2 CHARACTERIZING MULTIPLICITY AND SURVEY METHOD- 
OLOGY 

2.1 Key multiplicity properties 

The frequency of multiple systems and the companion frequency (average number of 
companion per target, which can exceed 100%), MF and CF respectively, are the most 
immediately quantifiable aspect of multiplicity. The difference between these two quantities 
hinges on the frequency of high-order multiple systems. Unfortunately, these two quantities 
are frequently confused in the literature, leading to ambiguous interpretations. Beyond 
these global quantities, and although the diversity of multiple systems is extreme, their 
statistical properties can be summarized in a handful of distributions, which we briefly 
introduce here. While correlations between these distributions could be present, multi- 
dimensional analyses require very large samples of binaries, which remain out of reach in 
most cases. 

The orbital period, P, of a binary system can be readily estimated for spectroscopic 
binaries (SBs). For visual binaries (VBs), the orbital semi- maj or axis, a, can be estim ated 
from the projected separation, p, on a statistical basis (e.g., Brandeker et al. 20061 ) . A 



power law parametrization, f(P) oc P a , is commonly used, w ith the part icular case a = — 1 



(a logarithmically-flat distribution known as "Opik's law", lOpikl 1 19241 ) being frequently 
favored. An alternative log-normal representation, parametrized by P and o"i og p, is also 
widely used, a lthough it differs only marginally from Opik's law for large values of ai os p 



( Heacoxl [19961 ). Distinguishing empirically between these two parametrizations would have 



fundamental implications on the binary formation process, as Opik's law suggests a scale- 
free process while a log-normal function implies a preferred spatial scale for companion 
formation. However, dynamical evolution throughout a system's lifetime can significantly 
alter its period, complicating the interpretation. 

The mass ratio, q = Af soc /Af pr i m < 1, can generally be derived from the observed flux 
ratio, although only a model-dependent estimate can be obtained for PMS and substel- 
lar objects. A power law distribution, f(q) oc q 1 , is frequently used to characterize the 
distribution of mass ratios and we use this formalism in this review to compare differ- 
ent populations. When no robust estimate of 7 has been proposed in the literature, we 
computed it from available data. The orbital eccentricity is an important tool to probe 
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the dynamical ev olution of a system: s hort period binaries can circularize as a result of 
tidal dissipation ( Koch fc Hrivnakl 198lh . while eccentricity can be pumped in three-body 
internal (Kozai) and external interactions. The expected distribu tion for a dynamically 
relaxed ("thermal") population is /(e) = 2e (|Ambartsumianl 119371 ) . Because eccentricity 
can only be reliably estimated for binaries whose orbit is short enough to be fully mapped, 
only a subset of all systems can be included in the analysis of the eccentricity distribution. 
We limit our discuss of this quantity to Section l5.1.4l Finally, the other orbital parame- 
ters define the orientation of the orbital plane relative to our viewing geometry and are 
therefore not as astrophysically compelling. The orientation of orbital planes is essen tially 



random with respect to our line of sight and to the stellar rotation spin vector (e.g.. lHale 



1994 ; |Popovic. Pavlovic fc Ninkovidl2004), while multiple syste ms only have a modest pref- 
erence for internal coplanarity (e.g.. Isterzik fc Tokovininl 2002[ ). We will not discuss these 
quantities further in this review. 



2.2 General methodology 

An optimal multiplicity survey should 1) rely on a complete, uniform search of a large, 
volume-limited sample and 2) make use of a variety of observing methods, as each of them 
is tailored to a certain type of companions. Magnitude-limited surveys, which are very 
common, are affected by an inevitable Malmquist-like bias, but this effect can be corrected 
statistically. Still, such surveys can be affected by more subtle selections biases, such as 
preferences for less extincted and/or more massive (brighter) targets. A sample size of at 
least 100 targets is typically required to measure multiplicity frequencies with a precision of 
±5%. Studies of orbital parameters (e.g., P or q) require a similar number of companions, 
hence a sample that is several times larger. 

Whenever possible, the comparison between two independent surveys should be per- 
formed directly over a c ommon section of the parameter space, and binomial statistics 
should be preferred (see Brandeker et al. 2006). However, most multiplicity surveys are 
incomplete and it is normal practice to estimate a completeness correction to take "missed 
companions" into account. Bayesian inference techniques can also be used to optimally 
exploit i ncomplete d ata while properly presenting the resulting degeneracies in parameter 
fits fsee lAlIer]|2007f ). but they require assumptions about parametric functional forms. If 



completeness corrections are taken into account, it is important to ensure that different 
surveys use very similar or identical assumptions. In addition, because of the statistical 
nature of incompleteness corrections, only robust estimates of CF can be obtained since it 
cannot be determined around which primaries are the companions missed in the survey. 



3 MULTIPLICITY ON THE MS 

In this section, we review the multiplicity of Population I MS stars grouped by mass bins: 
solar-type stars (Af* « 0.7-1.3 Mq; spectral types F through mid-K, Section ETTj) . low-mass 
stars (M* « 0.1-0.5 Mq; M0-M6, Sectionl33j> . VLM stars (M* <O.1M ; M8 and later, 
Section [3~3)) . intermediate-mass stars (M* « 1.5-5 Mq; B5-F2, Section [3~4)) and high-mass 
stars (M* > 8Mq; B2 and earlier, Section [33j) . The mass dependence of stellar multiplicity 
is addressed in Section loTTl 
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3.1 Solar- type stars 



3.1.1 Historical overview Solar-type stars represent an ideal population to probe 
stellar multiplicity on the MS. They are frequent enough to yield large, statistically robust 
samples and, while they are relatively bright and easy to observe, their luminosity contrast 
with any stellar companion is relatively modest. Solar-type stars were w ell represente d 
in the samples of bright nearby stars studied over the last century (e.g., Heintz 19691 ). 
but completeness and sample homogeneity were poor. Even the large sca le, uniform and 
dedicated solar- type multiplicity study conducted by Abt fe Levy 1 1976! ) had a limited 
completeness and suffered from the biases inherent to its magnitude-limited nature. 

The first modern, volume-lim ited multiplicity survey of solar-type stars was conducted 
by Duauennov fc Mavoi 1 199ll ). who studied 164 objects out to 22 pc by collating results 



from many separate studies and obtaining their own new radial velocity measurements. 
For the last two decades, this study has represented the "gold standard" of multiplicity 
surveys among field stars, against which most subsequent multiplicity surveys compared 
their results. Nonethe less, signifi cant incompleteness po tentially affected the conclusions of 
Duauennov fc Mayor , prompting Raghavan et al. 1 201(t ) to revisit this topic with a volume- 
limited sample of 454 star out to 25 pc. In addition to a much cleaner sample, that survey 
built on much improved observational methods to reach a very high completeness rate 
(« 97%). This study now represents the most complete multiplicity study of a well-defined 
sample and likely will remain a standard reference for years to come. 

In addition to these global surveys, which consider companions at all orbital periods, 
valuable results have been obtained in large single-method surveys which enable more uni- 
form detections rates. In the context of field solar-type stars, dedicated surveys for SBs, 



VBs and common proper motion compa nions have been conducted (e.g.. lHalbwachs et al 
20031 . iMason et al.lfl99sl iTokovininlboilh . 



3.1.2 Multiplicity frequency While Abt fc Lew ( 19761 ') concluded that "single 
stars are infrequent", both selection biases and survey incompleteness led these authors 
to significantly overestimate the compani on frequency of solar-type st ars. Instead of an av- 
erage of 1.4 companion per primary star, Duauennov fc Mavorl Jl99lf) and Raghavan et al 
( 2010h demonstrated that C^o^L 13M = 62 ± 3%. In addition. [Raghavan et all esti- 

44 ± 2%, i.e. 



mated that MF"j_ 13M@ = 44 ± 2%, i.e., a sl ight majority of all field solar-type stars 
are actually single, as is the Sun. Furthermore, iRaghavan et al.1 took advantage of their 



large sample size to identify a dependency of multiplicity on stellar mass. Specifically, they 
found that that super-solar dwarfs have a marginally higher multiplicity rate than sub-solar 



Similarly, their sample leads to 



dwarfs: MF^_\ 3 M@ = 50 ± 4% vs MFffi^ Mq = 41 ± 3%. 
J±^= 75±5% and CF n M 7 s _, „ = = 56 ± 4%. 
Mason et al.1 ( 19981 ) and iMetchev fc Hillenbrand! (2009j) have suggested that young (r 



CF? 



0.1-1 Gyr) solar- type dwarfs have an excess of visual companions over their older coun- 
terparts. However, these surveys were skewed towards super-solar targets, reducing the 
amplitude of the proposed excess. Only a bias-free survey could conclusively demonstrate 
that young field solar-type stars host more stellar (and substellar) companions than their 
older MS counterparts. 



3.1.3 Orbital period distribution The distribution of orbital periods for solar- 
type dwarfs has l ong been kno wn to be extremely broad, spanning up to 10 decades 
in orbital periods I Luvten 1930h . While Opik's law is a reasonable first approximation, 
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the distribution cuts off at both extremes (jHeintz 196S ), although the exact location of 



the "breaks" are somewh at ill-defined (e.g., Poveda. Allen fc Hernandez- Alcantara! 2007 



Tokovinin fc Lepinel 2012 ). Indeed, the high completeness and large sample size of modern 
volume-limited samples now allows to confidently exclude Opik's law as a good representa- 
tion of the orbital period distribution over the entire r ange of periods. On the other hand, 
the log- normal description that was first proposed by Kuiperl ( 1935b!) has n ow become a 
clearly better description of the observed distribution ( Raghavan et al. 20101 ) . The maxi- 
mum of the orbital period distribution occurs at P « 250 yr (a w 45 AU), with a dispersion 
of CTi g p ~ 2.3. 



3.1.4 MASS ratio distribution The difficulties associated with the detection of low- 
mass companions has lead to widely discrepant conclusions regarding the overall mass ratio 
distribution of solar-type binaries. In turn, it has been propos ed that this dis tribution is flat, 
bimodal, or monotonically increasing toward low-q systems I Trimble! 1990l . and references 
therein). Modern vol u me- li mited surveys have clarified this ambiguous picture. Although 
Duquennov fc Mayor! (| 1991 ) found a sing l e pea k around q ~ 0.3, the much higher com- 
pleteness of the survey by Raghavan et al. ( 20ld ) firmly established that this distribution is 
instead flat down to q « 0.1 with the exception of a marginally significant peak at q > 0.95. 
We fit a power law to the overall sample from iRaghavan et al] (adding all three panels in 
their Figure 16) and obtain 7^7-1.3 m @ = 0.28 ± 0.05. 

At a more detailed level, it is apparent that the mass ratio distr ibution differs significan tly 
between short- and long-period binaries. Splitting the sample of Raghavan et al. ( 2O10l ) at 
the median orbital period, short-period binaries are characterized by a strong peak at q ~ 1 
and a slowly declining f(q) function towards low mass ratios, while lon g-period binaries have 
a sing le peak around q fa 0.3, similar to the total distribution found bv lDuquennov fc Mayor 



( 19911 ). Specifically, we fin d that 7 n 7 



MS, log P<5.5 



(see also I Tokovinin! 12011 



= 1.16±0.16 and 7o 1 7 s i 1 ° S 3^> 5 - 5 = -0.01±0.03 



although neither power law is a perfect model of the observed 
distribution. 

One remarkable feature of the mass ratio distribution for solar-type binaries is the 
dearth of substellar companions, the so-called "BP desert ", first identified among SBs 



( Grether fc Lineweaver!l2006l : iMarcv. Cochran fc Mavorll2000l ) and further proposed for VBs 
(e.g.. iLeconte et ~ 2010| ). We discuss this in the general context of "extreme" mass ratio 
systems in Section [5^4l 



3.1.5 High-order multiple systems Triple and higher-order systems represent sa 
25% of all solar-type multiple systems, with a distribution of systems with n compo- 
nents that roughly follows a geometric distribution N(n) oc 2.5 _n up to n = 6 (see also 
Eggl cton fc Tokovin in 2008). However, this distribution is noticeably steeper for non-single 
systems: N(n) oc 3.7"". 

Not surprisingly, multiple systems are hierarchical. In all systems, the r atio of the orbital 
periods Piong/fshort has a mini mum value of > 5 (e.g.. lTokovininlll997l ). consistent with 



long term stability arguments (jEggleton fc Kiseleval Il995l ). The mass ratio distribution 
for the short period subsystems is flat up to q ~ 0.9, with a strong peak around q ?a 
1 which is largely res ponsible for the marginal excess reported in the overall sample of 
Raghavan et al.l ( 201C ). On the other hand, th e long period subsystems preferentially have 
a < 0.5 jRaghavan et aflboid iTokovininlkoOSh . 

More than half of all pairs with log P < 2 are members of high-order multiple systems 
and this proportion rises even further for the shortest-period systems (|Allen et al]|2012l . 
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Tokovinin et al]|2006h . In addition, the frequ ency of hierarchical triple systems amon g wide 
systems (separation >1000AU) is also high (jMakarov. Zacharias fc Hennessvll2008h . This 
suggests that the presence of an outer companion plays an important role in the formation 
of tight pairs, presumably through energy and angular momentum exchanges. It is plausible 
that this phenomenon is also related to the formation of equal-mass subsystems. 



3.2 Low-mass stars 



3.2.1 Historical overview The high incidence of low-mass stars enables large sam- 
ples of nearby objects to be searched for stellar companions. After th e first attempts a t 
detecting companions via the astrometric perturbation methods (e.g., Lippincott 197&t ). 



a series of instrumental breakthroughs over the last two decades have helped e stablish 
the mult iplicity of low-mass stars throug h radial velocity ( Marcv fe Benitz ll98S), direct 
imag ing ( Skrutskie. Forrest fc Shurelll989l ) and speckle interferometry ( Henry fc McCarthy 
1990) surveys. Although none of t hese surveys was v olum e-limited and their samples only 
partially overlaped, they allowed Fischer fc Marcv ( 19921 ) to perform a complete analy- 
sis of M dwarf m ultiplicity for 166 star s that has long been used as the low-mass stars 
counterpart to the iDuquennov fc Mavorl (|l991T ) survey for solar-type stars. More recently, 



similar-sized volume- limited surveys of low-mas s stars have been conducted (|Delfosse et al 



20041 . iDieterich et al.|[2012l . iReid fc Gizid[i997ah . In particular, the RECONS consortiun 
has gathered a quasi-complete sample of stars and companions out to ~ 10 pc (jHenrv et al 



200' 



3.2.2 Multiplicity frequency While most studies agree that low- mass stars have 
fewe r companions than solar-type s tars, published multipli city frequencies range from 26 ± 
9% jLeinert et al.lll997l ') to 42 ±9% ^Fischer fc Marcvll 1992V This spread can in part be ex- 
plained by the small number of companions detected in these surveys, further compounded 
by significant incompleteness corrections in certain domains. Modern volume-limited sur- 
veys, with their near-perfect completeness, draw a much more coherent picture, where 



MFi 



MS 



1 -O S M Q 



26 ±3% and CF T_ 



33 ± 5% l|Delfosse et alJbooi 



2012, Reid fc GizisHl997ar ). The widely-used higher frequency reported bv lFischer fc Marcv 



Dictcrich et al 



( 1992]) can be traced back to their analysis of common proper motion binaries, which focused 
on the earliest spectral types (M3.5 and earlier). Indeed, recent studies have demonstrated 
that the frequency of wide companions (50-10 5 AU) is a steep function of primary mass 
( Dhital et al.ll2010f ) so that an average over all spectral types leads to a much reduced com- 



panion frequency at large separations. Despite this trend, there is no significant variation 
of the total multiplicity frequency between early-M and late-M dwarfs since most of their 
companions are found at separations shor ter than 50 AU. O n the other hand, the frequency 
of spectr oscopic companions estimate d bv lFischer fc Marcvl has been upheld by more recent 
surveys (jdark. Blake fc KnaprJl2012r i 



3.2.3 Orbital period distribution iFischer fc Marcvl (|1992T ) suggested that a log 



normal distribution similar to that found for solar-type binaries could also apply to low-mass 
st ars. As pointed ou t above, however, the frequency of wide companions was overestimated 
bv lFischer fc Marcvl and the orbital period distribution must thus be revisited. A log-normal 
function is still acceptable for separation below 500 AU, with a ~ 5.3 AU and ai og p ~ 



1 http:/ /www. recons.org/ 
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1.3 based on the latest update to the RECONS sample. At larger s eparation, an Opik - 
like distribution seems to hold out to 0.2-0.3 pc for early-M primaries (jPhital et alJboid ). 
Considering the small number of companions included in volume-limited samples, and the 
difficulty to draw a bias-free picture from the other available surveys, it remains unclear 
whether and how these two distributions are physically connected. 



3.2.4 MASS ratio distribution The flat mass ratio distribution found bv lFischer fc Marcv 



( 19921 ) suffered from potential biases as only primaries in the M0-M3 range and companions 
significantly above the substellar limit were considered. Nonetheless, including lower mass 
M-type primaries and all secondaries yields a mass ratio distribution that is flat or slightly 
declining tow ards low-q syst e ms. F or instance, we have fit a power law distribution to the to- 



tal sample of lDelfosse et al.l (|2004l . their Figure 4) and found that 7 



1-0.5 Alp 



= 0.39 ±0.23. 



However, the situation is more subtle, since the mass ratio distribution is clearly depen- 
dent on the primary mass and the binary separation. Like solar-t ype binaries, short pe riod 
M-dwarf binaries are significantly biased towards high-g system s llReid fe Gizis 1997a ). In 



addition, low- 



ting the RECONS sample, we find that ■y^ ,a - 5 Au 



systems are limited to the earliest M-type dwarfs (|Dieterich et al.ll2012l ). Fit- 

0.5 Alp ~ 1 - 1 - x - u tlllu 70. 1-0.5 AIq — — 0.3±0.3. 
Similarly, we find that 7^3-0.5 a/ = -0.2 ± 0.3 and Yo*i-o.3M Q = L 9 ± 1-7. Finally, 
we note that BD companions are easier to find around M dwarfs than around solar-type 
stars. Although dedicated searched have discovered a number of such companions (e.g. 



Burningham et al.l 120091 : iForveille et"afl 1200-J : iNakaiima. Qppenheimer fc Kulkarnil Il995l : 



Radiean et al 



20081 ). the frequency of such companions is much lower than that of stel- 



lar companions based on volume-limited survey, with only 3 of the 23 companions in the 
RECONS sample being substellar. 



3.2.5 High-order mul tiple systems The ratio of multiple to binary systems among 
low-mass stars is ~ 21% (|Reid fc Hawlevl 120051 ). similar to the ratio observed for solar- 
type stars. Furthermore, the distribution of systems with n components among M dwarfs 
follows an N(n) oc 3.9~™ for n > 2, closely matching the distribution for solar-type sys- 
tem^. Finally, high-order multiple systems appear to have a much high er occurrence rat e 
among systems which contain the widest pairs among M-dwarf systems ( Law et al. 2O10l ). 
suggesting that higher system masses are needed to produce the widest systems. Overall, 
besides the lower overall multiplicity of M dwarfs compared to solar-type stars, the basic 
statistics of high-order multiple systems do not seem to differ significantly between these 
two populations. 



3.3 VLM stars and substellar objects 



3.3.1 Historical overview VLM objects, even the nearest ones, are typically very 
faint and therefore difficult targets for both high-resolution imaging and spectroscopic ob- 
servations. Broad molecular b ands and fa st rotational velocities further hamper searches 
for SBs. In his review of BDs, iBasril (|200ol ) only briefly addressed the issue of multiplicity 
as the first BD binary systems had just been discovered. In the following decade, a slew 
of high-resolution surveys identified ma ny systems, and thorou gh re views of multip licity in 
the substellar regime were presented bv lBurgasser et al.l 1 20071 ) and iLuhmanl (2012j). 



2 See for instance the VLM binary database (http://www.vlmbinaries.org/) 
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3.3.2 Multiplicity frequency In the absence of a large volume-limited survey for 
both VBs and SBs, a piece-wise approach is needed to estimate the total multiplicity 
frequenc y of VLM stars. The broad range of observed frequency of spectroscopic com- 
panions jBasri fc Reiners 2006 . Blake. Charbonneau fc White 20ld . Guenther fc Wuchterl 



20031 . iReid et alJbooj 



Tanner et al.l 



20121 ) is primarily a result of different completeness 



limits for each survey. Combining all these surveys, 77 objects of spectral type M8 or later 
have been observed at least three times with high-resolution spectroscopy, revealing 4 con- 
firmed SBs and at least two more candidate binaries. The frequency of SBs (a < 1 AU) 
is therefore at least 5.2jl^"f %. At larger separations (a > 2AU), high-resolution imaging 
surveys indicate an obse rved binary frequency of about 15±3% after correction for their 
inherent Malmquist bias (|Bouv et al ] |2003l . Burgasser ct al 1 120071). It is likely that a handful 
of companions remain undetected in the 1-2.5 AU range (e.g., Bernat et al. 2010|). Finally. 



no m ore than 2.3% of VLM objects have a companion in the 40-1000 AU range (jAllen et al 
20071 ) . Combining these various studies yields a total multiplicity fraction for field VLM 



objects of 2 0-25%, consistent with the frequency of CF^^q 1Mq 



22±4% estimated by 

Alien! (|2007h from a Bayesian analysis of several input surveys. We adopt this latter es- 
timate, which takes into account survey incompleten ess. Furthermore, while a few triple 



VLM systems have been confirmed or suspected (e.g.. lArtigau et al.ll201ll . iBurgasser et al 



2012), too few of these have been discovered to draw a meaningful conclusion. Instead, 
we simply estimate that CFf£~% x M ~ ^-^m»<o x M . While the studies discussed so 
far generally focus on late-M or L-type objects, recent surveys found a similar multiplicity 
frequency for the coldest (T-type) BDs. Even the higher frequency of visual companions 
among "transiti on" objects (late-L to early-T dwa rfs) can be explained by an enhanced 
Malmquist bias ( Burgasser. Looper fc Ravner 2010h . 



3.3.3 Orbital period distribution While too few VLM binaries are known to derive 
robust estimates, their orbital period distribution is much narrower than that of higher 
mas s objects. Most compa nions are found with projected separations in the 1-10 AU range 
(see Burgasser et al. 2007fh Only a handful of systems with projected separations larger 



than 50 AU are known i Billeres et al. 20051 . Burningham et al. 2010l . Dhital et al 



and th is near-absence of wide s yste ms cannot be accounted for by incompleteness (jAller 
2007T I. ~" 



Maxted fc Jeffries! j2005l ) and lAllenl (|2007l ) have used statistical analyses to derive 



(truncated) log-normal analytical descriptions of the overall orbital period distribution. 
Both studies agree regarding the peak of the separation distribution (a w 4-7 AU), but differ 
on the width of the distribution. It is probable that an intermediate width (ciogP ~ 0.4-0.5) 
is a good match to all observations. However, the orbital period distribution may follow a 
different analytical prescription and may even be asymmetric. More thorough searches for 
SBs are needed to settle this issue. 

3.3.4 MASS RATIO distribution The mass ratio distribu tion for VLM star jf| is heavily 
skewed toward equal mass systems (see IBurgasser et alj|2007i ) . The vast majority of known 



Biller et al 



20061 . 



systems have q > .7 and the most "extreme" systems have q ~ 0.5 (e.| 
Kasper et al.ll2007h . Several studies have attempted to reproduce the observed mass ratio 
distribution by assuming an intrinsic power law distribution and taking into account selec- 



* For substcllar objects, estimated masses depend on the assumed age of the objects and are therefore 
not as precise as for hydrogen-burning stars. However, the uncertainty applies similarly to both compo- 
nents of a presumed coeval binary, so that absolute uncertainties in q arc substantially smaller than on 
the masses of individual components. 
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Leggett 



tion and incompleteness effects (| Alien! 120071 burgasser et al.ll2006l . iLiu. Dupu-y 

While they all agree that 7a/„<o.i a/ q ^ can be firmly excluded, the estimated range 
of estimated values is very broad. This divers i ty is most likely the result of methodological 
differences. The shallow index found by I Allen! (|2007l ) is unlikely to be appropriate as it leads 
to the p rediction of many binar ies with q < 0.7, which is at odds with deep imaging surveys. 
Instead, Burgasser et al. ( 2006h argues that the sharp decline in the observed distribution 
of mass ratios as q decreases is only modestly amplified by the Malmquist bias. While this 



issue is not entirely settled, we adopt the steep distribution proposed bv lBurgasser et al 



3.4 Intermediate-mass stars 

3.4.1 Historical overview Intermediate-mass stars are much more challenging for 
multiplicity studies than lower-mass stars. The steep mass-luminosity relationship makes 
low-mass companions hard to detect, and given typical distances, there is an un-probed 
gap between SBs and VBs. Further complicating the picture is the presence, among nearby 
field A dwarfs of some ~ 30% of chemica lly peculiar stars. Some of these categories owe 
their peculiarity to multiplicity ( Abt 1965), inducing severe selection biases. In this review, 
we consider A stars as a unique population, averaging over the various sub-samples. 

Sea rches for SBs among i n termediate-mass stars hav e been ongoing for s everal decades 
(e.g., lAbt fc Snowdenl Il973l . ICarauillat fc Prieur! l2007l . ICarrier et afll2002l ). Systematic 
searches for visual companions have been largely incomplete for a long time, although 
some of the studies listed above have attempted to include wider companions. The efforts 
to identify low-mass visual companions have accelerated in recent year s as the presence of 
BP, and more rece ntly planetary-mass, companions is investigated (e.g. 
~ 20061 ). 



Ivanov et al 



Balega et al.ll2011 



In parallel to these studies of field stars, another avenue to probe the multiplicity of 
intermediate-mass stars is to consider rich stellar associations in the process of dissipating 
in the field. By an age of 5Myr, intermediate-mass stars already have reached the zero-age 
MS. The Scorpius-Centaurus OB association is the best stu died region to date, both for SBs 



Brown fc Verschueren 



2007al . IShatskv fc Tokovinin 



1997b and visual com panions (|Ko uwenho ven. Brown fc Kaperl 
120021 ). which allowed lKouwenhoveii"'et"luT " |2TO7b[ ) to perform 



a systematic analysis to derive the global multiplicity properties of intermediate-mass stars 
in that population. 



3.4.2 Multiplicity frequency The obs erved frequency of SBs among field intermediate- 
mass stars lies in the 30-45% range (|Abtlll983h . In the Sco-Cen OB association, the fraction 
of SBs is at least 30 %. The most thorough search for visual companions to A stars is the on- 
going VAST survey ( De Rosa et al.ll201ll ). In the separation range 50-2000 AU, over which 
the survey is complete down to the substellar limit, they find a frequency of companions 
of 40±4%, although many of these companions have q < 0.1 and some of them may be 
unrelated background star s (R. De Rosa, priv. comm.). In the Sco-Cen OB association, 
Kouwenhoven et al.l (|2005r ) found a frequency of visual companions of 37%. Considering 
the limitations of current surveys, the Sco-Cen population provides the most robust es- 
timate of the multiplicity of intermediate-mass stars. Current observations support the 
hypothesis that field stars have similar properties. Based on the various surveys listed 
above, we conclude that Mi ? j M 5 s _ 5 M > 50%, although incompleteness precludes robust 
estimates of the frequency of high-order multiplicity. In addition, the statistical analysis of 
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Kouwenhoven et all |2007bh resulted in CF^f-a m = 100 ± 10%. 



3.4.3 Orbital period distribution The orbital period distribution for field intermediate- 
mass m ultiple syst ems extends from less than 1 d to projected separations of well beyond 
10 4 AU (jAbtHl983h . Contrary to lower-mass field binaries, however, the orbital period dis- 
tribution for A stars appe ars to be bimodal, with a peak among SBs around P ~ 10 d (e.g., 
Carquillat fc Prieurll2007l ) and one for VBs around 350 AU (R. De Rosa, priv. comm.). In 
addition, there is marginal evidence for a dip ar ound P ~ 3 00 d, with a tertiary peak or 
a broad plateau for separations in th e 1-50 AU dBudailll999l). Even though the observed 
distribution seems more complicated, Kouwenhoven et al.l 1 2007b ) tried both a power law 
and a log-normal representation for the period distribution of the Sco-Cen intermediate- 
mass binaries with moderate success. In the former case, Opik's law is favored, while in 
the latter formalism, they find ai og p > 2.3. However, neither is a perfect description of the 
observed distribution and no simple analytical representation matches all observations. 



3.4.4 MASS ratio distribution The mass ratio distribution of field intermediate- 
mass stars remains elusive as a result of incompleteness and selection biases, which are 
largely responsibl e for the early claims of a strong preference for high-q systems (e.g., 
Abt fc Levvll 19851 ). Taking those limitations into account, it now appe ars that the mass ra- 
tio d istribution for SBs is either relatively shallow (e.g., 7 = — 0.3±0.2: ICarquilIat fc Prieur 
2007 ) or can be described by a broad Gaussia n distribution centered around q ~ 0.4 
(jCarrier et all I2OO2I . IVuissoz fc Debernardil 120041 ). Among visual companions, the VAST 
survey suggests a mass ratio distribution 7 ~ —0.6 that is roughly consistent with that 
found for SBs (R. De Ro sa, priv. comm.). Fin a lly, w e not e that a similar mass ratio di stri- 
bution was proposed bv lShatskv fc Tokovinin (|2002l) and iKouwenhoven et al.l (|2007bl ) for 



UpSc 



= -0.45 ± 0.15. Since this 



the Sco-Cen association. The latter study derived fi 1 5 _ 5MQ 
analysis includes both SBs and VBs, we adopt it as representative of all intermediate-mass 
multiple systems. 



3.5 High-mass stars 

3.5.1 Historical ove rview The multiplicity of high-mass stars has long been the 
subject of intense scrutiny ( Zinnecker fc Yorkej2007l . and references therein). However, their 
large distances from the Sun, large rotational velocity and the extremely high brigh tness 
contrast requirements conspire to draw an incomplete picture IjSana fc Evans! 120111 ). In 
addition, most high-mass stars are found in their birth environment, massive clusters or 
OB associations, and only ~ 20% of all nearby O stars are found in the field, many of w hich 
are thought to be ejected "runaway" stars (e.g., Chini et al. 20121 . de Wit et al. 2004). As 
a result, studying only "field stars" is likely to yield a very biased picture; 



Hig h-mass SBs have been studied in all environments for several decades (|Abt. Gomez fc Lew 



1990l : IChini et al.ll2012l : iGarmanv. Conti fc Massevlll980h . Individual clusters and OB as 



sociations, despite smalle r sample sizes, are more appropriate for a unifor m characteriza- 
tion of multiple systems (jKiminki et alj 12012 : ISana. Cosset fc Evans! 120091 ). Furthermore, 
it has recently become possible to probe the multiplicit y of high-mass sta rs in the low- 



metallicity environment of the Large Magellanic Cloud ( Sana et al. 2012tJ ). Surveys for 



visual companions to high-mass stars have been conducted both for "all-sky " samples (e.g. 
Mason et al.ll20o"9.lTurner et al. 20081 ) and single clusters/associations (e.g., Duchene et al 



2001 



Peter et al. 



I boiaY 
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3.5.2 Multiplicity frequency The remarkably broad range of estimated frequency 
of spectroscopic companions among high-mass stars is due to a combin a tion of small sample 
sizes, selection biases and incompleteness. Most recently, I Sana et al.l (|2012af ) estimated a 
bias- and completeness-corrected frequency of spectroscopic companions of 70±9% out to 
P ~ 3000 d and down to q ~ 0.1. Early-B stars have a lower SB frequen cy than their 
higher-mass counterparts (52 ± 4% for B2-B6 primaries; Chini et al. 20121 '). Due to the 
diversity of wavelengths and observing techniques, imaging surveys yield an even broader 
range of frequencies. A consensus has emerged that the f requency of visual companions 
over two decades in projected separations is ~ 45 ±5% fe.g. JPeter etaflboij iTurner et al.l 
20081 ) . Incompleteness in these surveys is modest down to q ~ 0.1 overall but very uncer- 
tain below this limit. There is sugge stive evidence that the frequency of vis ual companions 
decreases with the mass of primary (|Peter et alJboil IPreibisch et al.|[l999l ) . Interpolating 
between the populations of SBs and VBs suggests a " missed" comp a nion frequency in the 
10-20% range, in agreement with the small survey by iNelan et al.l (120041 ). This leads to 



9>0.1 



130 ± 20% and CF, 



9^0.1 



100 ± 20%. The latter is substantially lower 



C 1 F' 

I " " ' O 

than the estimate from lAbt. Gomez fc Levvl (|1990f ) but the large incompleteness correc- 



tion of that survey probably explains this difference. In addition, the number of high-mass 
stars that are truly si ngle stars is remarkably low. Considering only confirmed companions, 



MF 



M*>16 Mq 



80% (|Chini et al]|2012l , lMason et al.ll2009l . lSana fc Ev"ans]|201lh. The situa- 



tion is less well constrained for early-B stars, and MFs-i&m q > 60% (|Abt. Gomez fc Lew 
1990) is a conservative lower limit since many tight visual companions were missed in that 
survey. High-mass stars in rich clusters and OB associations have essentially identical mul- 
tiplici ty frequencies, which in turn is sign i ficantly higher tha n that of field and runaway 
stars jChini et al.ll2012l , iMason et aljlioog . ITurner et al.ll2008h . 



3.5.3 Orbital period distribution The relative lack of SBs with published orbital 
solutions has long impaired definitive analyses of the orbital period distribution for high- 
mass multiple systems. The high frequency of short-perio d binaries proposed early on 
( Abt. Gomez fc Levvl [l99(il ; iGarmanv. Conti fc MasseyHl980P ) was potentially subject to a 



strong selection bias, but has b een confirmed by modern surveys, which find a pile-up of bi- 
naries with 4 d < P < 8d fe.g.. lSana et al . 2012a). Power law fits suggest a slowly-declining 
(a « —0.5) distribu tion out to P = 3 000 d for both Galactic and Large Magellanic Cloud 
high-mass binaries ( Sana et al. 2012al lbh . However, extrapolating even a moderately steep 
distribution results in far too few visual companions, suggesting a more complex functional 
form. One possibility to match all observations combines two independent distributions: 
a population of short period binaries (logP < 1, 30% of all high-mass stars) and a power 
law period distribution extending out to > 10 4 AU. While an Opik-like law for the sec - 
ond population is compatible with available observations (|Kiminki fc Kobulnickvl 120121 ) . 
imaging surveys which repeatedly find a "peak" at th e shortest probed angular separations 
rather support a slowly declining d istribution instead ( Mason et al. 20091 . Nelan et al. 2004 . 
Peter et alJboil ISana et al.ll201lh . 



3.5.4 MASS ratio distribution Limited sensitivity to even intermediate-mass com- 
panions has long plagued studies of the mass ratio distribution of high-mass stars. Further- 
more, a large fraction of their c ompanions are still in the PM S p hase, resulting in m odel- 



dependent mass estimates. Both lKiminki fc Kobulnickv (2012) and lSana et al.l (12012 



eluded that the intrinsic distrib ution of mass ratio s for high-mass SBs is essentially flat (e.g., 
= — 0.1±0.6 for Sa na et al ] |2012al . which we adopt here). While both studies 



P<3000 d,g>0.1 
Tm*>16 M a 
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find a peak around q ~ 0.8, they rule out the presence of a se parate population of equal-mas s 
(so-called "t win") system s among SMC eclipsing binaries ( Pinsonneault fc Stanek 20061 ) . 
Among VBs, iPeter et al, |2012l ) concluded that a flat distribution and a broad Gaussian 
distribution centered on q ~ 0.45 match equally well their observations of high-mass stars 
in the Carina region . On the other hand, surveys of the Orion Nebula Cluster (ONC; 
Preibisch et alj|l999h and NGC6611 (jDuchene et al.ll200lf ) found a strong preference for 



low-mass companions (q < 0.5), although some of these systems may not be physically 
bound. Combining these three samples and fitting a power law to the resulting distri- 
bution yields J^F^wMq ~ ~0-55 ± 0.13, indeed favoring lower-mass companions but not 
significantly different from the distribution found for SBs down to ~ lMg. Lower mass 
companions remain mostly out of reach of present surveys but the upcoming generation of 
high-contrast, high-resolution imaging instruments will soon provide a first foray into this 
category of "extreme" systems (q « 0.1). 



4 MULTIPLICITY OF OLD AND YOUNG STARS 

In this section, we consider multiplicity properties of stellar populations of old MS stars 
(Population II stars) and at increasing earlier stages of stellar evolution, from open clusters 
to the youngest embedded protostars. The time evolution of multiplicity is addressed in 
Sectionl5T2l 



4.1 Population II stars 

The multiplicity properties of Population II stars can be used to probe the star formation 
process in low metallicity environments. While early studies led t o conflicting conclusions 
regarding a possible deficit of comp anions to Population II stars ( Abt fc Willmarthl 1987 . 



regarding a possible aencit 01 companions to copulation 11 stars l|ADt &c winmartnlliaofl . 
Jaschek fe Jaschek 19591 . Port 19261 ). a more coherent picture has recently emerged thanks 



to instrumental improvements and larg er and better-defin ed samples of low-metallicity 
dwarfs and high-velocity halo stars (e.g.. ICarnev et al.lll994T l. 

SBs among Population II stars appear to be equall y frequent and to share very similar 



prop erties with their higher-metallicity counterparts (jGoldberg et al.l |20o3 . lLatham et al 



2002). On the other hand, the frequency of VBs (a > 10 AU) is markedly lower than for 



solar -metallicity stars ( Allen. Poveda fc H crrcra 2000: Lodieu et al.ll2009l ; Zapatero Osorio fc Martm 



2004), though perhaps not for K-M dwarfs (|Jao et al.ll2009l ). As a result, the total multiplic- 
i ty frequency of Population II sta rs is CFq j^_ \ 7 3 M = 39 ± 3 % and CFq_°^ 6 m = 26±6% 
( Jao et al.ll2009l . lRastegaevll2010l . respectively) . 

The distribution of orbital period among solar-mass Population II multiple systems is 



characterized by a 
period tail out to 



narrow peak in orbital period around log P_ 



2-3, accompanied by a long - 



10" 4 AU that follows roughly Opik's law (| Allen, Poveda fc Herreralboool) 
and a marked dearth of systems with orbital periods shorter than lo g P ~ 1 ([Latham et al 



2002 j RastegaevlboiCll ) . The secon dary peak at ~ 50 AU proposed by Zinnecker. Kohler fc Jahreifj 



(|2004r ) has not been confirmed bv lRastegaev (2010 , casting doubt on the its existence. Too 
few low-mass Population II binaries are known to robustly estimate their orbital period 
distribution, though they may represent a broader distribution than for solar-metallicity 
stars (Jao et al. 2009). Finally, while the mass ratio distribution of solar- type Population II 
binaries is still poorly determine d, companion masses down to the substellar limit are found 
in roughly uniform distribution ( Zapatero Osorio fc Martml 2004 ). 



Gaspard Duchene & Adam Kraus 



4.2 Open clusters 

Members of open clusters represent young (typically 50 Myr to 1 Gyr) and uniform pop- 
ulations that will eventually dissipate in the galactic field. These clusters are dynamically 
old, in the sense that their ag e is much larger than their crossi ng time, and wide binaries 
have long been disrupted (e.g., lKroupa|[l995l . barker et al]l2009l ) . Despite the advantage of 
large and homogeneous samples, multiplicity studies of open clusters are surprisingly scat- 
tered in the literature, with little overlap between samples. The best studied open clusters, 
in which both SBs and close VBs have been searched for, are the nearest ones (D < 200 pc): 
a Per, the Pleiades, Praesepe and the Hyades. No dramatic difference has been identified 
between these clusters in terms of their multiplicity properties so that a composite average 
can be constructed to draw a global picture. 



4.2.1 Multiplicity frequency We first address the multiplicity of solar-type stars, 
which are the best-studied in open clusters. Most su rvey for SBs in open clusters agree on a 
13-20% frequency o f spectroscopic companions (e.g.,|Griffm 2012. IMermilliod. Grenon fc Mayor 
20081 . |Prosseilll992h . Although adaptive optics surveys rely on substantial incompleteness 
correction for the closest com panions (p < 1"), th ey are more sensitive to low- mass com- 
panions than speckle surveys (jMason et al.l Il993al lbf. batience et~aH Il99sl . I2002T I and their 
companion frequency of ~ 25% over 1.5-2 decades in separation is more robust. 

Combining SB and VB surveys therefore yields an observed companion frequency of open 
clusters of ~ 40% for solar-type stars, although some companions most lik ely still escape 



detect ion. Taking advantage of the range of distances to various clusters, iPatience et al 



( 2002h found a total companion frequency of 48 ± 5% for all separations up to 580 AU 
and with mass ratios q > 0.25. Assuming a mass ratio distribution similar to that of field 
FGK dwarfs (7 ~ 0), the total frequency is therefore CF° p 7 e *f^ e * ~ 65%. Independently, 
photometri c studies have suggested multiplicity fr equencies in the 65-70% range for solar- 
type stars (|Converse fc Stahler|[20oj . lKahlerlll999h . 

Low-mass cluster members have only been systematically studied via spectroscopic an d 
imaging surveys m the nearby Hyades cluster (jReid fc Gizislll997bl . lReid fc Mahonevll200ol ). 
Combining these two surveys leads to CFq i^s"!/™ ~ 35 ± 5%. Finall y, VLM star s and 
BDs in open clust ers have only been searched for visual companions (|bouv et all bood 



Martin et all 120031 ). Extrapolating the observed fre quency of visual com panions assum- 
ing an overall Gaussian orbital period distribution, Gizis fc Reidl (1995) suggested that 
q F^an c uster — 2T ± 16%. Photometric studies, which may be affected by a mass seg- 
regation selection effec t , broadly support these results despite some significant scatte r 
jBoudreault et alJboiAlConverse fc Stahlerlliooj . lLodieu et allbooi IPinfield et al.ll2003l ). 



4.2.2 Orbital period distribution Although no single cluster has been studied 
over the entire range of orbital periods, it is possible to estimate the average orbital period 
distribution in open clusters in a piece-wise manner using the various surveys discussed 
above. T he distribution for so lar-type stars is broad and unimodal over the range < 
logP < 7 ^Patience et alj|2002l ). The second peak at logP < 1 proposed bv lGriffinl (|2012l ) 
likely stems from the fact that their sample include s intermediate-mass st ars. 

The median separation of ~ 4 AU estimated by IPatience et al.l (|2002l ) based on a log- 
normal fit, somewhat tighter than among field stars, may be a consequenc e of incomplete ness 

among VBs. Indeed, the fr equency of VBs included in their analysis by Patience et al.l ap- 

pears to be underestimated ( Bouvier. Rigaut fc Nadeaull997 ; Bouvier et al Jl200ll :l Morzinski ct al 
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2012), suggesting that the median separation is significantly larger. It is therefore probable 
that the orbital period distribution of solar-type binaries in open clusters is indistinguish- 
able from that of field stars. Unfortunately, chance projection of unrelated members and 
background stars prevent any analysis of systems wider than a ~ 1000 AU in open clusters. 



4.2.3 MASS ratio distribution The distribution of mass ratios among solar-type bi 



narie s in open clusters is essentiall y flat for bot h SBs (e.g., Bender fc Simonl2008l . Goldberg fc Mazeh 



1994 . iMermilliod fc Mavorl 1 19991 ) and VBs (|Patience et alj I2002T ). In a ddition, we note 



that substellar companions t o solar-type stars in open clusters are rare (lltoh et al. 2011 
Metchev fc Hillenbrandll2009l ), although a few have been recently discovered IjGeifller et al 



20121 . iRodnguez et al.ll20120. On the other hand, multiple systems are increasingl y skewed 



towa r d high-q for low-mass jReid fc Gizis|[l997bl . lReid fc Mahonevl2000h and VLM |Bouv et al 
20061 . lLodieu et aT]|2007l . iMartm et al.ll2003l ) members of open clusters. 



4.3 Young nearby associations 

The past decade has seen the discovery of several "moving groups", i.e., groups of co- 
moving, coeval young stars that are unb ound and not associated with a recognizable site 
of star formation (see Torres et al. 2008t ). The ages of these groups (~ 77-100 Myr) make 
them a valuable bridge between the sites of ongoing star formation on one hand and open 
clusters and the field population on the other. Most of these groups are small (N < 50), 
however, so they individually do not provide robust statistical constraints on the full bi- 
nary population. Furthermore, many of the known binaries were discovered serendipitously 
rather than in statistically robust surveys. The few analyses to date suggest that they are 
broadly consistent with what is seen in the field and in yo ung star-forming regions (e.g., 
Brandeker. Javawardhana fc Naiitall2003l ; Evans et alj|2012l ) , excep t for an anomalously low 



frequency of wide binaries (a > 20 AU) in the r\ Cha moving group I Brandeker et al. 20061 ) . 
Because of their youth and close proximity to the Sun (d < 50 pc), mem bers of these associa- 
tions have been subject to surveys aiming primarily at planet discovery ( Chauvin et aljBoiol . 
and references therein). These have revealed a small number of planeta ry-mass companions 



l Chauvin ct al 



Lowrance et al 



2004 



1999 



Lagrange et al] 200E ) and numerous BDs (e.g., (Bowler et all I2OI2I , 



Neuhauser fc Guent hcr 2004. 1Wahhai et alJl20Tlf ). as w ell as provid- 



ing r obust constraints on the existence and extent of the purported BD desert (| Evans et al 
20121 ). A uniform analysis of all these results is still missing. 



4.4 PMS stars 

We now turn to multiplicity surveys targeting populations of PMS stars with ages in 
the 1-5 Myr range and with stellar mass M* < 2 Mq , collectively called T Tauri stars 



(TTS ), or Class II/III sources depending on whether they host a protoplanetary disk (|Lada 



19871 ) . Multiplicity surveys have targeted TTS located in a broad range of environments: 
dense stellar clusters (such as the ONC and IC348), young OB associations (e.g., Sco-Cen), 
and "loose" (so-called T) associations (like the Taurus-Auriga or Chamaeleon star-forming 
regions). The latter populations typically contain many less stars (~ 100) than clusters 
and OB associations and only a handful of intermediate-mass stars. The closest of these 
populations are located ~ 150 pc away from the Sun. 
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4.4.1 Historical Overview Even though the existence of close pairs of TTS was 
noted very early with seeing-limited observations ( Herbid 1962 . Jov fc van Briesbroeck 
19441 ). the multiplicity of TTS was only addressed by dedicated surveys in the 1990s. 
Since then, thorough multi- wavelength surveys of nearby star-forming region have allowed 
the definition of large, complete samples for multiplicity studies (e.g., Herbig fc Bell 1 19881 , 



LuhmarJIioTj ). In addition to seeing-limited surveys |Reipurth fc Zinneckerl 1 19931 ) and 



mining of all-sky databases (jKraus fc Hillenbrand! l2009al ). closer companions have been 



searched with a variety of high-angular resolution techniq ues, down to 2-3 AU resolution. 



In m o re or less chronolog ical order, speckle int erferometry ( Ghez. Neugebauer fc Matthews 



1993: iLeinert et al.lll993l). lunar occult ation |Simon et al 



19951). Hubble Space Telescope 



imaging ( Padgett. Strom fc Ghe zll 19971). adaptive optics I Duchene. Bouvier fc Simonlll999l ) 



and aperture masking I Kraus et al.l 200a ) have all been used in this context. The search for 
tighter companions remains challenging as strong and variable emissi on lines hamper radial 
veloc ity analyses. Nonetheless, a small number of SBs are known ( e.g., Mathieu. Walter fc Myers] 
19891 ). The only systematic surveys published to date |Meloll2003l , lNguven et al.ll2012l ) found 
a companion frequency similar to that of field stars but yielded no new orbital solutions. 
In the remainder of this section, we only discuss VB surveys, which have been much more 
extensively studied. 



4.4.2 Multiplicity Frequency Binary surveys of PMS stars in T associations gen- 
erally have measured com panion frequen cies that are roughly twice as high as among 

all 



eraily nave measured companion trcqucncics tnat are rougniy twice as mgn as ami 
solar- type MS stars (e.g., IDuchene 19991 ). This led to the widespread meme that 



young stars are born as multiples. Th e largest complete surve ys to date have been con- 
ducted in Taurus and Upper Scorpius ( Kraus et al. 20081 . 2011). In Taurus, for instance, 



CF, 



Class II/III,Tau 



— 64ig 1 % and Cfo^J^Mg 11 " = 791^% over the separation range 3- 
5000 AU. Despite these high frequencies, it appears that ~ 1/4- 1/3 of all sy s tems consist 



of a single star, revealing the prevalence of high-order multiples. iKraus et alj (|201ll ) noted 



that the distribution of degrees of multiplicity was broadly consistent with a purely Poisso- 
nian distribution, N(n) oc e~ n , and hence that adding additional companions resembles a 
stochastic process. This distribution resembles the geometric distribution seen for the field 
(Section 3.1.5; Eggleton & Tokovinin 2008). 

At the other extreme, numerous studies have shown that the multiplicity frequency in 



young dense clusters is much lower. This applies to the locally extreme ONC ( Kohler et al 



2006 . iPetr et al.lll998T l. but also to other young cl usters in the Orion complex jBeck. Simon fc Close! 
2003 ) as well as to IC 348 in the Perseus cloud ( Luhman. McLeod fc Gol densonll 2005l). In 
the ONC, where the largest sample of targets has been assembled, Rerpurth et al.l ( 2007 ) 
found ci?ciassii/m,ONC = 8.8±1.1% over a decade in projected separation, a factor of ~ 2.5 
lower than found in Taurus over the same range. Although it has long been proposed that 
multiplicity is a smooth function of the stellar density, the statistical precision of current 
surveys in in termediate regions, such as the Ophiuchu s cloud, do not conclusively support 
his scenario jKing et alj|2012al . iMaxks fc Kroupall2012l ). 



Although they lon g remained unacce ssible, several surveys have targeted BDs and VLM 
stars in recent years (jBiller et al.ll201ll . and references therein). Collectively, they demon- 
strated that most young low-mass objects have no companions at projected separations 
>5-10 AU. Furthermore, Kraus fc Hillenbrand |2012l ) performed a combined analysis of all 
published surveys and showed that the binary frequency declines precipitously from 0.5 Mq 



to below the substellar boundary, with CF 



Class II/III 
Af+<0.1 M Q 



< 11% (50% confidence) 
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4.4.3 Orbital Period Distribution As for the field population, both a log-normal 
and power law have been used with reasonable success to describe the distribution of or- 
bital periods among TTS for the last two decades. The inability to distinguish between 
these distributions results in part from the limited range probed by observations. Recent 
studies suggest a separation distribution that follows Opik' s law out of ~ 500 AU for solar- 
type stars and slowly declining for lower-mass stars fe.g.. iKraus et alj|2011f ). As a result, 
the global distribu tion is consistent with both a very bro ad log-normal, wit h a\ og p > 2 
(jKraus et al.l 120111 ). or a slowly declin ing power law (e.g.. King et al.ll2012bT l. The latter 
conclusion also matches the results of Kouwenhoven et al. ( 2007bl) for intermediate-mass 
multiple systems in Upper Scorpius. One important caveat is that t he binary population 
in dense clusters drops sharply at pr ojected separations of >200 AU (jReipurth et al.ll2007l ; 
Scallv. Clarke fc McCaughrearj[l999h . 

The orbital period distribution of binary TTS appears to be truncated, wit h very few sys- 
tems outside of a maximum separation that varies smoothly with stellar mass (jKraus fc Hillenbrand 
201 21 . and references therein) . The truncation radius appears to be ~ 300 AU for stars in the 
0.25-0.5 Mq and as small as ~ 25 AU for the lowest mass TTS. However, the paucity is not 



complete: a small numb er of "unusually wide" pairs have been identified (|Konopackv et al 



20071 . IKraus et al.ll201ll . and references therein). These are discussed in a broader context 
in Section [5~5l 

Finally, we note that the distribution of orbital separations is skewed towards tighter 
systems in disk-less Class III objects relative to Class II objects which host circumstellar 
disks. This is likely a result of disk evolution, since tidal forces in tight bi naries (a < 50 AU 



can severely truncate disks and result in much shorter disk lifetime ( Cieza et al. 200 



can severely trun 
Kraus et al]l2012h 



9, 



4.4.4 Mass Ratio Distribution Given the shallow mass-luminosity relation for PMS 
stars, imaging surveys are quite sensitive to low-mass companions, routinely reaching q ~ 0.1 
and mo re recently probing down to q < 0.01. However, uncertainties in PMS evolutionary 
models (jHillenbrand fc Whitdl2004h and the possibility of circumstellar dust emission and 
differential extinction could work to smooth out and/or bias the genuine mass ratio distri- 
butio ns, particularly for photometrically estimated mass ratios (jWoitas. Leinert fc Kohlei 
2001) . Nonetheless, in the absence of comple te spectroscopic follow-up campaigns (e.g., 

cases. 



Hartigan fc Kenvorj|2003l . Iwhite fc GhelbOQll ) , this is the only available tool in most i 

Most sur veys of solar-type binary TTS have found m ass ratio distributions that are 
roughly flat (jKraus et al ]|2008l . 2011, Lafrcnicrc et al jgoosh . Using a power law parametriza- 



tion for the mass ratio distribution, Kraus et al.1 1 201lh found To'^-^s™. — 0-2 ± 0.2 for 



the entire Taurus population, albeit with a marginally significant trend with stellar mass, 
whereby lower mass systems preferen tially have q > 0.5. This tren d is further expanded by 



considering the lowest mass binaries: Kraus fc Hillenbrand ( 20121 ) found T^o^^i^Af ~ 1- 



Due to the paucity of truly substellar binary systems, it is uncertain whether this trend 
continues through the substellar regime. We note, however, that there is a small fraction of 
objects that seem to host binary companions that are much less massiv e than the primary, 
and may also have unusually wide orbital radii (e.g., Luhman et al.|[2~009l ). These are further 



discussed in Section[5~5l At the other extreme, intermediate-mass binaries have a small but 



significant preference for low-q systems (|Kouwenhoven et al.ll2007br ) 
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4.5 Embedded protostars 

Over the last decade, increasing attention has been devoted to studying th e multiplicity of 



embed ded protostars, objects in the earliest phases of stellar evolution (e.g.. IDuchene et al 
2007b! ) . These represent the best chance of probing the formation mechanism of multiple 



systems at its origin, before subsequent internal and external dynamical evolution can affect 
unstable systems. However, protostars also pose serious challenges to companion searches. 
By definition, they are still embedded in the remnants of their nascent core and hence are 
best studied at infrared and radio wavelengths, where high-resolution techniques became 
available only recently. Furthermore, the mass of the central objects are extremely difficult 
to estimate. 

Protostars are traditionally categorized based on their spectral energy distribution. Class 
sources are the least evolved objects, with more mass in circumstellar material, mostly in the 
form of a dense envelope, than in the central object (jAndre. Ward-Thompson fc Barsonv 
19931 ). In contrast, the envelope mass of Class I sources is lower than that of the central 



object. We discuss both categories separately below. 



4.5.1 Class I protostars Although they are deeply embedded, Class I sources are 
bright near- and mid-infrared s ources. To date, about 350 Class I source s have been observed 
with seeing-limited imaging ([Connellev. Reipurth fc Tokunagal I200S L 



Haisch et al.1 



20041) 



and, for about 100 of them, with high-resolution imaging techniques (|Connellev. Reipurth fc Tokunaga 
20091 . IDuchene et"Zll2007al ) . Until a global statistical analysis of these various surveys is per- 



formed, we estimate the multiplicity of Class I in the following manner. First of all, we only 
consider sources within 500 pc of the Sun and companio ns with flux ratios AL < 4 mag to 
limit i ncompleteness issues. The direct imaging survey of lConnellev. Reipurth fe Tokunaga 
( 2008|) yields 27 com panions in the 200-2000 AU range out of 136 targets, while the adap- 
tive optics surveys of Duchene et al. ( 2007a ) and Connellev. Reipurth fc Tokunaga ( 20091 ) 
found 15 companions between 50 and 200 AU to 88 targets (with additional companions 
at shorter separations for the closest targets). Combining these two fractions, we find a 
total frequency of companions to Class I sources over the 50-2000 AU separation range of 
CF GlassI = 36.9 ± 5.2%. The associated distribution of projected separations is consistent 
with Opik's law out to 5000 AU although c hance projection of unrelated protostar s could 
be non- negligible at the largest separations ( Connellev. Reipurth fc Tokunaga 20081 ) . 

There is marginal evidence for an environment dependency among Class I protobinaries 
in that the distribution of separations for Orion s hows a marked peak in the 60-200 AU 
range not seen in othe r nearby star-forming regions (jConnellev. Reipurth fc Tokunagall2009l . 
IDuchene et al.|[2007al ). In this context, it is important to emphasize that Class I sources are 
distributed through the various molecular clouds in the Orion complex (a majority are in the 
L1641 cloud) and none belongs to the ONC. In other words, even the Orion Class I sources 
studied to date were born in dynamically quiet environments, so that any difference in that 
population is more likely to be a consequence of some other environmental difference than 
dynamical evolution. There is also evidence for a temporal evolution of the properties of 
protobinaries, as the most embedded, and presumably youngest, systems almost exclusively 
possess wide companions (1000 AU an d wider), the frequency of which rapidly declines as 
one considers more evolved protostars (jConnellev. Reipurth fc Tokunagall2008l ). 

The search for much tighter companions among Class I protostars is remarkably chal- 
lenging as radial velocity measurements are hamper ed by extremely weak pho tospheric 
features. Based on 2 or 3 measurements per systems (jViana Almeida et al.ll2012l ) or com- 
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paring a protostar radial velocity to that of gas in its parent cloud ( Covev et al, 20061 ) . 
small-scaled surveys suggest that some systems tighter than a few AU have already formed 
by the Class I phase, but their interpretation remains ambiguous. 



4.5.2 Class PROTOSTARS The rarity of Class sources makes it necessary to group 
all observed targets in a single sample to achieve meaningful statistics, precluding studies 
of environmental dependencies. In addition, high-resolution observations at nea r-infrared 



and c entimeter wavelength s are probing scattered light and free-free emission (jReipurth 



200Cl . iReipurth et al.ll2004f ). respectively, neither of which can be used to unambiguously 
detect companions. Instead, the optimal observational probe of Class multiplicity is there- 
fore millimeter (broadly speaking, 850 /j.m-7 mm) interferometric mapping, even though 
the linear resolution typically achieved with such observations is only ~ 100 AU. Further- 
more, clustering on scales of 1000 AU to 0.1 pc prevents studies of the widest systems (e.g., 
Carrasco-Gonzalez et al]|2012l ; iPeretto. Hennebelle fe Andral2007l ). 



Early multiplicity studies of Class protostars suggested that they have a very high mul- 
tiplicity freq uency although their target se l ection was mos tly haphazard and potentially 
biased (e.g.. lLoonev. Mundv fc Welch! lioool ; iReipurthlboOCh . Recent high-resolution mil- 
limeter surveys, wh ich primarily focus ed on the disk/envelope structure, reported low com- 
panion frequencies ( Maury et al.l 20 ipt). albeit with extrem ely small sample size. Adding to 
the sample surveyed by M^uTvetalT Enoch et al. 1 20 111 ) concluded that only 1 out of 11 
targets had a companion within 2000 AU (albeit with minimum separations probed ranging 
from 50 to 400 AU). However, a much larger number of Class sources have been imaged at 
high spatial resolution, leading to the discovery of several other c ompanions in that separa- 
tion range, thereby casting doubts on their apparent rarity (e.g., Loonev. Mundv fc Welch! 
200Cl ). Further complicating the picture is the possibility that some continuum point 



sources could be related to o utflow features instead of being bona fide companions (e.g., 
Maury. Qhashi fc Andrej|2012T) . a possibility that can only be clarified with follow-up obser- 
vations. A global re-analysis of all high-resolution millimeter maps of Class sources in the 
context of multiplicity is crucially missing but beyond the scope of this review. For now, 
we simply note that the multiplicity of Class sources is not statistically significantly lower 
than that of the more evolved Class I protostars in the ~ 100-2000 AU range. 



5 OVERALL TRENDS AND "EXTREME" SYSTEMS 

We now gather the results discussed in the previous two sections to highlight the main 
trends associated with multiplicity as a function of primary mass and age. We then discuss 
in more detail several remarkable categories of companions: twin companions, extreme 
mass ratio systems and extremely wide systems. Possible physical implications of empirical 
trends are discussed in Section[6] 



5.1 Multiplicity as a function of stellar mass 

In this section, we focus on the mass-dependence of multiplicity properties on the MS 
(see Section|3] and Table[T}. With continuously improving completeness and larger sample 
sizes, our empirical knowledge will soon allow decisive tests of the predictions from star 
formation theories. 
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Table 1: Multiplicity properties for Population! MS stars and field BDs 



Mass Range 


Mult./Comp. 


Mass Ratio 


Orbital Period 




Frequency 


Distribution 


Distribution 


AL < 0.1 M 


MF = 22t\ % 
CF = 22 + _\ % 


7 = 4.2 ± 1.0 


Unimodal (log-normal?) 
a fa 4.5 AU, criog P fa 0.5 


0.1 M© < M* < 0.5 M Q 


MF = 26 ± 3 % 
CF = 33 ± 5 % 


7 = 0.4 ±0.2 


Unimodal (log-normal?) 
a fa 5.3 AU, criogp fa 1.3 


0.7 M < M* < 1.3 M Q 


MF = 44 ± 2 % 
CF = 62 ± 3 % 


7 = 0.3 ±0.1 


Unimodal (log-normal) 
a fa 45 AU, criogp « 2.3 


1.5 M < M* < 5M Q 


MF > 50% 
CF = 100 ± 10 % 


7 = -0.5 ±0.2 


Bimodal 
Pfa lOd & o fa 350 AU 


8 M < M^ < 16 M Q 


MF > 60% 
CF = 100 ± 20 % 






M* > 16 M© 


MF > 80% 
CF = 130 ± 20 % 


7 p<30ood = _o.i ±0 .6 

7 a>100AU = _ 05±01 


Peak + power law 
P fa 5d & Opik (?) 



5.1.1 Multiplicity frequency As has long been known, the overall multiplicity 
frequency of MS stars is a steep, monotonic function of stellar mass (see FigureQJ. Given 
the limited number of points available and the fact that companions to high-mass stars with 
q < 0.1 are not accounted for, we refrain from proposing a functional form that reproduces 
the intrinsic dependency of CF(M*). The dependency of MF(M*) is qualitatively similar 
to that of CF(M*) except for a shallower slope resulting primarily from the fact that MF 
has an intrinsic upper bound at 100%. 

The frequency and order of multiple systems also are steep functions of ste llar mass. For 

insta nce, all 11 sextuple systems in the Multiple Star Catalog are A-type stars (lEggleton fc Tokovinin 



20Q8), while a single such candidate system is known among nearby solar- type stars (jRaghavan et al.l 



2010j) , and the highest-o rder multiple system containing only low-mass stars currently 



known has 4 components (|Reid fc Hawlevll2005h . While it appears that the ratio of binary 
to higher-order systems does not vary much for objects with M* < 1.5 Mq, it probably rises 
significantly toward high-mass stars. 



5.1.2 Orbital period distribution The distribution of orbital periods is unimodal 
for solar-type and lower-mass stars, but both the median separation and width of the 
distribution decrease sharply with decreasing stellar mass. As a result, the frequency of 
companions in the 1-10 AU range (10-15%) does not vary significantly with stellar mass 
for M* < 1.5 Mq, including substellar objects. The much lower overall multiplicity of field 
BDs traces a marked deficit of wide companions as opposed to a uniform depletion over all 
separations. 

Intermediate- and high-mass stars have more complex distributions of orbital periods. 
A strong peak at the shortest periods (logP fa 0-1) is found in both populations with 
an amplitude that increases with stellar mass. VBs show a peak for intermediate-mass 
stars, while the situation is less clear for high-mass stars, for which a shallow power-law is 
currently preferred. Interestingly, the frequency of companions in the 1-10 AU range among 
intermediate-mass appears in reasonable agreement with that observed among solar-type 
stars. 
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Figure 1: 

Dependency of CF (red squares) and MF (blue triangles) with primary mass for MS stars and 
field VLM objects. The horizontal errorbars represent the approximate mass range for each 
population. For B and O stars, only companions down to q as 0.1 are included. The frequencies 
plotted here are the best-estimate numbers from Sections 3.1-3.5, also reported in Tablcfl] 



5.1.3 MASS ratio distribution Although a simple power law representation is imper- 
fect for most samples, this formalism offers the most straightforward criterion to compare 
multiple systems of various masses. As shown in FigureO the observed distribution of mass 
ratios is close to a flat distribution ( | < 0.5) down to q ~ 0.1 for all ma sses M* > 0.3 Mq, 
extending to high-mass stars the conclusions of lReggiani fc Meven (|201ll ). Below this limit, 
the mass ratio distribution becomes increasingly skewed towards high-g systems. Further- 
more, shorter-period systems among solar-type and low-mass stars have a steeper mass 
ratio distribution than wider systems. 

Many past studies have favorably compared the mass distribution of companions to that 
expected from random pairing from the IMF for field objects. In this situation, a nd leaving 
aside the "smearing" induced by broad ranges of primary masses (|Toutlll991f ). 7 should 
increase from -2.3 (Salpeter's slope) at the high-mass end to ~ around or below the 
substellar limit, which does not match observations at any primary mass. Although a 
proper comparison between observations and theory /simulations should be conducted on 
the raw f(q) distributions rather than on estimated power law indices, the hypothesis of 
random pairing from the IMF is robustly excluded by observations (see Figure[2)). On the 
other hand, the rarity of substellar companions to low-mass stars suggests an alternative 
model in which the mass ratio distribution is flat between q — 1 and a minimum companion 



that is independent of the primary mass (at least up to M* ~ 1.5 Mq) and 
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Figure 2: 

Power law index fitted to the observed distribution of mass ratios for multiple systems as a 
function of primary mass. Red diamonds represent fits to the overall population of multiple 
systems within a certain range of primary masses, whereas blue squares and green triangles 
represent fits to the subsets of "tight" (P < P) and "wide" (P > P) binaries, respectively. 
Horizontal crrorbars represent the mass range for each subsample. The dashed cur ve represe nts 
the in dex that would be derived if the companions followed the single stars IMF of IChabrierl 
l|2003l l and a simple power law would be fit to the resulting companion mass distribution over the 
mass ratio r ange 0.1 <<?<!■ The a sterisk marks the derived power law index using the single 
stars IMF of lBochanski et alj l |2010h . The power law indices plotted here are the best-estimate 
numbers from Sections 3.1-3.5, also reported in Tablc[T] 



lies around or somewhat below the substellar limit. However, given current statistical 
uncertainties, it is also possible that there is no such minimum companion mass, or that it 
varies with primary mass. 



5.1.4 Eccentricity distribution Early studies of the eccentricity distribution i n 



stellar binaries were hampered by severe selection biases (e.g.. lAitkenlll932l . iHeintz 196S ) 



More modern systematic surveys have now yielded a clear and uniform picture (e.g., Abt. Gomez fc Levy 



199C| : lDe Rosa et al-lboiilPupuv fc Liul201ll;|Kimmki fc Kobulnickvl2012l ; lKonopackv et al 



20ld ; lRaghavan et alj|20ld ; ISana et alj|2012al ). 

The eccentricity distribution for field multiple syst ems shows r emark ably little depen- 
dency on primary mass, extending the conclusions of Abt (2005, 200d ) across the entire 
range of primary masses. Beyond a period of ~ 100 d, all populations show an essentially flat 
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Figure 3: 



Cumulative distribution of eccentrici ties for systems with 2 < log P < 4 for field multiple systems 
among solar-type stars (green curve; [ptaghavan et al.| [201(f). low-mass st ars (orange curve; fr om 
the SB9 catalog, iPourbaix et al.ll2004l) , VLM stars an d BDs (red curve; iDupuv fe Liull201lh , 
high-mass stars (blue curve; from the SB9 catalog and lAbjlioO^JSana et al.ll2012af) . The 
dot-dashed curves indicate incomplete samples, for which the eccentricity distribution is 
potentially biased. The dashed and dotted curves represent the expected distribution for a flat 
and thermal distribution, /(e) = 2e, respectively. 

distributiorff). As illustrated in Figure[3] a broad Gaussian distribution centered on e ~ 0.4 
is a possible alternative representatio n, since it produces somewhat less high-eccentricity 



systems (e.g.. IStepinski fc Bl ack 2001). Discriminating between these two representations 
hinges on large, unbiased samples, which are still missing for most stellar masses. In any 
case, the observed distributions are all inconsistent with the so-called "thermal distribu- 
tion" , even though the latter is frequently assumed in numerical simulations. 

Almost all short-period systems have circular orbits as a result of tidal dissipation. The 
observed circularization periods, which range from < 1 d to « 20 d, is a function of the 
primary mass and a monotonic function of stellar age, from open clusters members to 
Population II stars. A handful of exceptions with eccentric short-period binaries are known 
among solar-type binaries. All of them are either members of higher-order systems that may 
be u ndergoing Kozai cycles, and /or relatively young so that they have not fully circularized 
yet jMathieu fc MazehlHiij ). 



Between these two regimes, the eccentricity distribution has a well-defined rising upper 



4 Although Duqucnnoy & Mayor ( 1991) favored such a distribution, it clearly violates their bias- 
corrected distribution, which is fiat beyond e ~ 0.3 (see their Fig. 6b). 
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envelope in the P-e pla ne. The mild correlation between mean eccentricity and orbital 
period, first proposed by Finsen ( 19361 ) and borne out by recent surveys, is mostly due to 
the dearth of systems with e > 0.6 at orbital periods 10-300 day. Intriguingly, extrasolar 
planets populate the P — e plane in a very similar manner as stellar multiple systems 
( Udrv fc Santo3l2007l ). 



5.2 Multiplicity as a function of age 

We now turn to the temporal evolution of multiplicity for solar-type and lower-mass stars, 
which we can probe for ages ranging from < 0.1 Myr for the most embedded protostars to 
> lOGyr for Population II field stars. Higher-mass stars are only probed over a much 
narrower age range and do not allow for as thorough an analysis. 

5.2.1 Multiplicity frequency No population other than nearby field stars has been 
sufficiently studied to draw a complete picture of multiplicity, so this analysis can only be 
achieved for certain categories of companions. For SBs, no significant difference is observed 
between TTS, open cluster members, Population I and II field stars, although this conclusion 
is potentially weakened by incompleteness and biases and deserves further attention. For 
the remainder of this section, "multiplicity frequency" specifically refers to the frequency 
of visual companions. 

Even among VBs, only a narrow range of separations (~ 50-500 AU) has been studied 
across all populations, due to the diversity of distances involved and observational tech- 
niques employed. This limited range is insufficient to draw statistically robust conclusions, 
given typical sample sizes. Instead, we proceed by estimating for each population the com- 
panion frequency per decade of projected separation using the widest possible separation 
range within the approximate 10-3000 AU range. This implicitly assume that all population 
of multiple systems follow Opik's law. Although an imperfect fit to the overall observed dis- 
tribution, it is a reasonable approximation in the separation range considered here (which 
is near the peak of most proposed log-normal separation distributions) . 

The resulting frequency of visual companions is shown in FigurelH where we have elected 
to focus on the largest samples studied with uniform and complete techniques whenever 
possible. In addition, in cases where several surveys address the same "population" (e.g., 
several open clusters), we average the resulting multiplicity frequency and consider the 
dispersion of individual results as an additional source of uncertainty to account for some 
inherent diversity. Among TTS, we group separately dense clusters from loose associations, 
since their multiplicity frequency is markedly different. 

Comparisons between populations of young stars and that of field stars is complicated 
by the broad range of stellar masses probed in the early stages of stellar evolution, where 
surveys typically include objects of all masses from the substellar limit up to ~ 2 Mq . In 
the past, it has been customary to compare most PMS multiplicity surveys to the well- 
studied population of solar-type field stars. However, a fairer comparison should instead be 
performed with a weighted average of solar-type and low-mass stars that takes into account 
possible selection biases in the young stellar populations. 

Figured] highlights two parallel trends. On the one hand, stellar multiplicity does not 
vary on timescales < 50 Myr in low-density environment (star-forming regions and young 
associations). On the other hand, there is no significant difference in multiplicity frequency 
between young and open clusters on one hand, and Population I field star on the other hand. 
It remains undetermined whether the dichotomy of multiplicity properties at early ages nec- 
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Figure 4: 

Dependency of the frequency of visual companions per decade of projected separation with age for 
solar-type (~ 0.7-1.5Mq, blue triangles), low-mass (~ 0.1-0.5 Mq, red squares) and overall 
populations of young stars (~ 0.1—2 M p,, orange diamonds). Results from individual surveys are 
indicated as gray symbols and include [Duchene et al.l j2007a|) and 

Reip urth fc Tokunagal j200Sft for Class I so urces. iRatzka. Kohler fc Leinertl d2005h 



G 



oimcllcv. 



[ 200 71). iBeck. Simon 



, and lLafreniere et al. I <200ct> for T/OB associa tions. iReipurth et al,l 
Close! ll2003ll.lpuchene. Bouvier fc Simonl dl999h and 



lLuhman ] McLeoTT*fc"cio*ktensqnri 20051) for clusters of TTS, iBrandeker, Jayawar dhana fc Naiital 
l|2003l 'l: lBrandeker et alJ <2006ft. IChauvin et alj fcOlOT) andlMcCarthv fc White! <2012ft for nearby 
young associations. and lBouvier. Rigaut^T^adeaul jl99 7| ) nB^uvier*et*Tl* for open clusters. 

The horizontal errorbars represent the approximate age range for each population while vertical 
errorbars compound statistical uncertainties and the observed dispersion between similar regions. 
The dashed lines represent constant frequencies describing qualitatively the predicted behavior of 
low-density associations (top line) and dense clusters (bottom line). The hashed region represent 
the range of possible behaviors for stellar clusters in the embedded phases. 



essarily implies that the initial multiplicity properties are environment-dependent. This is 
discussed in more detail in Section[6] Finally, the marked multiplicity deficit among Pop- 
ulation II stars, which is limited to wide binaries, suggests either a metallicity dependence 
in binary formation or evolution on a timescale of several Gyr. 

5.2.2 Other multiplicity properties Beyond the evolution of the multiplicity fre- 
quency, most of the trends identified among PMS multiple systems closely match those 
observed among open cluster members and field stars. This includes the dependence of the 
mass ratio distribution on both stellar mass and separation, and the dependence of the max- 



Gaspard Duchene & Adam Kraus 



imum separation with stellar mass (e.g., Kraus et al. 201ll ). Possibly the most discrepant 
property when comparing young and mature multiple systems is the orbital period distri- 
bution: among young stellar objects it appears closer to Opik's law than among field stars 



espec ially for solar- type systems (e.g., Connellev. Reipurth fc Tokunaga 20081 : 



Kraus et al 



20 111 ). However, this only applies for young multiple systems in loose associations. Given 
the marked deficit of wide companions in young clusters like the ONC, averaging the orbital 
period distribution over all star-forming regions (weighted by their relative star formation 
rates) would result in a distribution that is much more similar to that observed among field 
stars. Not surprisingly, it thus appears as if multiplicity properties of field stars can be 
explained from those of TTS without much evolution. 



5.3 "Twin" binaries 

The existence of a distinct population of "twin" binary systems has been repeatedly 
advocated and is still subject to intense debate, in part because of the loose definition of this 
phenomenon. Broadly speaking, twin binaries are near equal-mass, short- period binaries 
which are over-represented over the underlying population of binaries (e.g., iLucv fc Ricco 
19791 ). 

An excess of systems with orbital period <43d (med i an period ~ 7d) and q > 



among 0.5-2Mq binaries has been detected (Lucy 20061 . Simon fc ObbielBooj ) although 
this subpopulation represents only 2-3% of all SBs. This low frequen cy explains why it is 
not readily apparent in volume-limited surveys ( Raghavan et al. 2010l ). It is probable that 
this excess results from accretion in the early stellar evolution stages that acts to equalize 

t he components' masses. 

Pinsonneault fc Sta nck (2006) found an excess of systems with q > 0.8 7 among ecli psing 



binaries among 7-27 Mq stars in the Small Magellanic Cloud. However, Lucv ( 20061 ) sug- 
gested that this trend could be explained in large part by selection biases and subsequent 
studies of Galactic high-mass stars do not find evidence for a separate population of twin 
systems (see Section [5^3)) . Still, it is plausible that twin binaries may be a feature asso- 
ciated with low metallicity, since high-mass binaries in the Large Magellanic Cloud show 



a mar ginally steeper mass ratio distribution than their Galactic counterpart (|Sana et al 
2012b! ) . 



5.4 "Extreme" mass ratio systems 

The interest in "extreme" mass ratio systems (broadly defined as q < 0.1) has rapidly 
grown, both in the context of the search for planetary-mass co mpanions a nd because of the 
difficulty for binary formation model to produce them (e.g., Bate! 20121 ). Recent surveys 
have largely overcome the contrast limitations of earlier searches, providing new insights 
into the formation of VLM companions. Figure[5] which is based on large-scale surveys 
of field and young stars, illustrates the marked deficit of low-g systems, especially among 
(very) low-mass objects. 



The extremely low frequency of BP companions among solar-type SBs (jGrether fc Lineweaver 
20061 : iMarcv. Cochran fc Mavojboool ) is inconsiste nt with extrapolations of both the binary 
mass function and the planet ary mass function (e.g. jHoward et al ]|20ld iReffert fc Quirrenbachl 
201 ll . ISahlmann et al] 120111 ). This suggests a sharp disco ntinuity in companion masses 



around the hydrogen burning limit ( Thies fc Kroupal 20081 ). although this BD "desert" is 
not completely arid and warrants further analysis. 
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Figure 5: 

Le/t; Distribution of mass ratios as function of mass for nearby field objects with JW* < 1.5 Mq. 
The dotted and dot-dashed curves indicate constant companions masses at the substellar and 
planetary regime limits (0.075 and 0.013 Mq), respectively. Data for solar-type stars (green 
diamonds), low- mass s tars (orange triangles) and VLM objects (red pl us signs) are from 
IRaghavan et al.l i20ld 'l. the RECONS survey and ljanson et al.l j2012al 'l. and the VLM binary 
database, respectively. BD companions to low-mass stars identified outside of lar ge-scale surveys 
are sh own as open orange triangles. The few systems with M pr ; m < 0.5 Mq from IRaghavan et alj 
correspond to lower mass subsystems in hierarchical multiple systems. Right: Similar plot 
for star-forming regions and young associations. Shown as gray diamonds in this plot are 
companions to Class II/IH objects in T aurus, Upper Scorpius and Chamaeleon I, using the l arge 
surveys from lKraus et al.l j2008l . 1201 J) , iKraus fc Hillenbrand! i2012l ), lLafreniere et all j2008f l. In 
addition, red open diamonds represent a number of objects discovered via small-scale surveys or 
pointed observations with a particular emphasis on VLM primaries and/or companions. 



Among visual companions, the stringe nt limits on the presence of substellar compan- 



ions obtained a decade ago (<1% ; e.g.. Iffinz etafl I2OO2I , [McCarthy fc Zuckermanl |20o4 



Oppenheimer et al.ll200ll ) applied to samples largely dominated by low- mass stars. More 



recently, deeper surveys focusing on solar-type stars found potentially higher fr equencies, al- 
thou g h the very small number of detections limits the robustness of this result JCarson et al 



200fj . iLeconte et aUboid iMetchev fc Hillenbrand |2009| , iTanner. Gelino fc Lawlboioh . It 



remains unclear whether these substellar companions are a simple extension of the p opu- 
lation of companions below 0.075 Mq, as proposed by IMetchev fc Hillenbrand or 
whether a distinct trough in companion mass is present at large separations as for short- 
period systems. One approach to test whether there is a demarcating line between two 
distinct populations of visual companions consists in studying the frequency of (very) low- 
mass stellar companions to intermediate- and high-mass stars, which probe similar mass 
ratios. In this context, w e note the existence of a handful of pla netary-mass objects around 



intermediate-mass stars (Kalas et al. 200. 



Marois et al 



200ctl but an apparen t dearth of 



companions around the substellar limit (jjanson et al l" l201ll . Ivigan et al.l l2012l ). Up 
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ing surveys for wide planetary mass companions to nearby stars using "extreme" adaptive 
optics instruments will help address this issue in a statistically sound manner. 

Multiple surveys have addressed the existence of the B P desert among visual companions 

in both star-forming region s and nearby associations (e.g., Evans et al.ll2012 ; Kouwenhoven. Brown 



2007al : lTanner et alj|2007h. These survey generally found that BP companions ar e rare but 



Kaper 



not absent (e.g., Janson et al. 2012bl . Lafreniere et al. 2011 . Nielsen et al. 20121 ). Finally, 
high-resolution imaging surveys also have discovered a small number of planetary-mass com- 
panions (~ 4-20 Mjuy) at distances of > 100AU from young stars in nearby star-forming 
regio ns ( Chauvin et al.11200 4 L Ireland et al Iboill lLafreniere. Javawardhana fc van Kerkwiikl 
20081 . iLuhman et alJbOQrj . lNeuhauser et al.f2005l . lSchmidt et al.ll2008h . Unfortunately, most 



of these discoveries were serendipitous, precluding a robust statistical analysis of their fre- 
quency. 



5.5 Extremely wide systems 



Observationally, the correlation between the median separation and primary mass ob- 
served for stars with M* < 2 Mp> is accompa nied by a similar trend for the maximum 



binary separation a max (e.g.. iReid et al.ll200lT ) 



an unconfirmed discontinuity at ~ 0.2 Mq have been proposed (Burgass cr et al 



Although differen t functional forms and 

200l . 



Close et al. 20031 . Reid et al. 2001), this upper limit on binary separation appears to hold 



for both field populations and P MS systems fe.g.. iFahertv et al.ll201ll . iKonopackv et al, 
2001 iKraus fc Hillenbrandl[2009al ) . Furthermore, several authors used energetic arguments 
to propos e that the maxim um binary separation is actually set as a minimum binding en- 
ergy (e.g. lciose et al.ll2003l ). It remains undetermined which quantity is the most physically 
relevant. 

A key to understanding the physical meaning of the a ma x — Mtot correlation lies in "unusu- 
ally wide" systems, which violate the general trend. While the vast majority of field binaries 
fall inside the maximum separation limit, numerous studies have successfully searched for 
systems well beyond that limit. For solar-type and intermediate -mass stars star s , sys- 
tems with separation out to several pc have been identified (e.g., Caballero 20091 . 2010l . 
Shava fc Ollinelboilh . Among field low-mass stars, for which the term "u nusually wide" 



was f irst coined, systems at separations > 10 3 AU have been found (e.g 



Artigau et al 



2001 ICaballero![2007l . bhital et al~ll2010l . IFahertv et allboill . iRadigan et al.ll2009l ). A 



ber of unusually wide low-mass systems have also been found in nearby star-forming r egions 



( Allers et al.ll2009l . lBeiar et all2008ljKraus fc Hillenbrandl2009al . lLuhman et alj|2009l ). Around 
the substellar limit, although a marked paucity of companions is observed beyond >25 AU, 
a small number of wide pairs have been identified, particularly among young objects ( e.g., 



Beiar et al . 2008, C hauvin et a l. 2004, Kraus 



HillenbrandlboollLuhman et al-ibood ). In 

all cases, these "unusually wide" systems have separations a > 10a max . Even though 
follow-up observations of some of these wide syste ms indicate that they are frequently part 
of hierarchical multiple systems ( Law et al. 2010l ). their binding energy is still well below 
the empirical limit discussed above. 

The existence of these unusually wide systems, coupled with the fact that the empiri- 
cal limit on binary separation has been repeatedly revised upwards over the last decade, 
suggests that there is no absolute upper bound on binary separation, but rather a gradual 
rarity beyond a certain separation. 
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6 IMPLICATIONS FOR BINARY FORMATION AND EVOLUTION 



Multiplicity is a smooth function of stellar mass, both on the MS and in the earlier 
phases of stellar evolution. This trend points to a similar formation process across a broad 
range of prestellar core masses. Core ("prompt") fragmentation still appears as the leading 
candidate scenario to form multiple systems, although the extreme diversity of systems, 
most notably the larg e breadth of t he orbital period distribution, may require at least two 
separate mechanisms 1 Tohline 20021 ) . 

Fragmentation appears to be mild, in that most cores only produce between 1 and 3 
components. Since the frequency of quadruple and higher-order systems among field stars 
is low, most high-order systems would need to be destroyed early on if core fragmentation 
typically led to a larger number of stellar seeds. This destruction would result in a much 
lower overall multiplicity frequency among field stars, especially among low-mass stars, 
since many ejected components would be single stars. While this is s upported by purely 
dynamical star formation simulations (e.g., Delgado-Donate et al.ll2004l ), it is not borne out 
by observations (see also lGoodwin fe Kroupall2005l \ Numerical simulations including the 
effects of magnetic field and/or radiative feedback appear to match more close ly observed 
multiplicity properties since these processes severely hinder fragmentation (e.g., Bate 20121 . 
Hennebelle fc Tevssierl [20081 ). 

The st atement th at "most stellar systems formed in the Galaxy are likely single and not 
binary" jLadalbood ) is still undecided. While it is true that most field stars are single as a 
result of the preponderance of low-mass stars, it is contrary to observations in low-density 
star-forming regions. It is unclear whether these regions are representative of multiple star 
formation in all environments. Indeed, the universality of multiplicity properties among 
different birth environments is still in doubt. There are marked observational differences, 
most notably the lower multiplicity observed in clusters relative to loose associations of 
young stars. However, numerical simulations sh ow that (internal and c luster) dynamical 



Marks fc Kroupall2012h . suggesting that 



evolution is by and large over after 1 Myr (e.g. 
current observations do not probe multiplicity in an early enough evolutionary stage to be 
conclusive. Because frequent interactions in dense clusters can reduce the initial multiplicity 
frequency by a factor of 2-4, a picture in which multiplicity f requency is universally high 



at early times remains plausible (jKroupa fc Petr-GozensH2011I ). One critical test would be 



to determine the multiplicity properties of clustered protostars since Class I surveys have 
so far only targeted objects born in "quiet", low stellar density environments. Another 
possible test of this scenario consists in comparing the mass-dependency of multiplicity in 
young clusters and T/OB associations, since dynamical simulations pred ict that lower mass 
systems are more prone to disruption (jKroupa. Aarseth fc Hurlevll200ll ). 

Dynamical evolution can naturally explain both the "low multiplicity" and "high multi- 
plicity" tracks observed in FigureSl in which clusters and loose associations represent ex - 
treme cases of highly effective and ineffective binary disruption 1 Kroupa fc Bouvier 20031 ) . 
Irrespective of whether initial multiplicity properties are universal, the existence of these 
two tracks has immediate implications on the origin of field stars. Indeed, it appears that 
most stars form in clusters, an argument already made in the past based on multiplicity 
( Marks fc Kroup"all2011 . Parker et al.ll200^ . Patience et al.|[2002h as well as from star/cluster 



counts ( Lada fc Ladaf 20031 ) . A more thorough analysis that goes beyond the multiplicity 
frequency and also considers other parameters, such as the orbital period distribution or 
dependencies on the primary mass, is still needed. 

The fact that the mass ratio distribution is essentially flat for most primary masses 
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(A/* > 0.3 Mq) but becomes increasingly steeper at lower masses fSection l5.1,3p sets im- 
portant constraints on both binary formation itself and the subsequent accretion through 
which final stellar masses are built. In particular, the suggestion that there may be a lower 
limit to companion masses, which remains unconfirmed with current observations, could 
be pointing to the smallest scale associated with fragmentation. In addition, the (small 
but growing) population of planetary-mass objects at large distances from their parent star 
(Section l5.4p poses a significant challenge to existing models of star and planet formation. 
Their orbital r adii are too large for a "plane tary" formation model (i.e., from a proto- 
planetary disk; IBqss! [200 ll . IPollack et al1ll996h . but their masses are extremely l ow for a 



'binary" formation model (fa lling near or below the opacity-limited minimum massi lKratter 



Murrav-Clav fc Youdin2010)r ). Nonetheless, it appears that one of these unlikely outcomes 



must occur. Determining whether these objects represent an extension of (sub)stellar com- 
panions to extremely low masses or form a separate population of objects will be addressed 
in the next few years. 

At short orbital periods, although the long-discussed BD desert has been repeatedly 
confirmed as tracing physically distinct formation mechanisms, the remarkable similitude 
between the eccentricity distributions of planetary and stellar-like companions (including 
massive BDs) is intriguing. One possible interpretation is that the final eccentricity in both 
planetary and multiple stellar systems is ultimately set by dynamical JV-body interactions. 
If this is true, it remains to be understood why the eccentricity distribution of multiple 
systems does not follow the thermal distribution, however. 

The observed correlation between primary mass and "maximum" binary separation is 



plausibly linked to the typical prestellar core size (both on the order of ~ 0.1 pc Bahcall fc Soneiral 



1981 . Lada et "al1l2008h . Similarly, the median binary separation, which also scales with pri 



mary mass, could be hinting at a preferred spatial scale for fragmentation. Alternatively, the 
physical origin of a max could lie in the dynamical evolution of systems. Indeed, the widest 
systems have a very weak binding energy that makes th em vulnerable to disruption by ot her 



nearby passing stars, either in their birth environment llKroupall 19951. 



Parker et al 



2009h 



m the Galactic field ( Jiang fc Tremainl 2010l : Weinberg. Shapiro fc Wasserman 19871 ). The 
latter secular evolution may be responsible for the lower multiplicity observed among Popu- 
lation II stars at large separations. A constant minimum binding energy can be expected in 
dynamically relaxed populations, such as virialized stellar clusters or field stars. However, 
even dynamically primordial environments like Taurus-Auriga and Upper Scorpius broadly 
follow the same a max trend, casting doubt on this interpretation. 

The physical origin of unusually wide systems (Section l5.5[l is not fully ascertained. They 
could form in a similar manner as other binaries from a single core. Since they would not 
survive dynamical interactions in dense clusters, one possible explanation for their extreme 
paucity in the field population is that "unusually wide" binaries form via a low-probability 
channel, and that only those b inary systems born in unbound associations can remain 
bound until they reach the field ( Kraus et al. 20081 . Kraus fc Hillenbrand! 2009a ) . The main 
Mtot — amax relation would then indicate the typical constraint from protostellar core sizes 
or the maximum spatial scale for fragmentation. Alternatively, they may have formed from 
separate cores and have been ejected s imultaneously from a cluster on sufficiently simila r 
trajectories as to remain weakly bound ( Kouwenhoven et al.ll201ol . Moeckel fc Clarkell2011 ). 
The low probability of this scenario readily accounts for their rarity among field stars. 



Stellar Multiplicity 29 



7 SUMMARY POINTS 



• The frequency and key properties associated with stellar multiplicity vary smoothly 
with primary mass, suggesting a similar formation process across all core masses. 

• Multiplicity properties appear to be mostly set at the PMS phase. Subsequent evo- 
lution is modest at best, as expected from dynamical evolution models. 

• No single population of binary systems has a mass-ratio distribution that agrees with 
random pairing from the IMF. Also, no population has an eccentricity distribution 
that follows the "thermal" distribution, /(e) = 2e. 

• Observations of both PMS and field stars do not support the existence of a firm 
upper limit on binary separation, but an increasing rarefaction towards the largest 
scales (e.g., ~ 1 pc for solar-type systems). 

• At short orbital periods (log P < 3), there is a marked distinction between planetary- 
and stellar-type systems, presumably tracing different formation mechanisms. 

8 FUTURE ISSUES 

• Are the multiplicity properties resulting from the star formation process universal or 
do they depend on the native environment? The mass-dependency of multiplicity 
among PMS in different environments stars holds a critical clue. 

• Are there significant differences in the multiplicity properties of high-mass stars rela- 
tive to lower-mass stars that would not follow a smooth mass dependency? A signifi- 
cant frequency of stellar companions with q < 0.1, or the confirmation of a separate 
but substantial population of "twin" binaries, would set high-mass stars apart. 

• What is the formation path and associated timescale for SBs? The frequency and 
properties of VBs appear to be set by the TTS phase, but is this also true for SBs? 

• What is the origin and fate of "unusually wide" systems? 

• Is the dearth of "extreme mass ratio" systems (known as the BD desert among tight 
solar-type and lower-mass binaries) independent of separation and primary mass? 
What are the physical processes responsible for this phenomenon? 

• Does multiplicity hinder or promote planet formation? How does it depend on binary 
separation, mass ratio, and eccentricity? 
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10 SIDEBAR 

Note to editor: place alongside Section 2 

Short period binaries (log P < 4) are best identified through radial velocity monitoring. 
Longer period ("visual") systems are revealed by direct and increasingly higher-angular 
resolution imaging surveys using speckle interferometry, adaptive optics-aided imaging, the 
Hubble Space Telescope, aperture masking on monolithic telescopes, lunar occultation and 
long-baseline interferometry. Companions in the most challenging "intermediate" regime 
(P > 5yr, p < 0.05") can also be probed by astrometric monitoring, although it has been 
difficult to apply to large, well-defined stellar populations except for bright stars observed 
by Hipparcos. 

Multiplicity can also be probed indirectly. In uniform stellar populations, high mass ratio 
systems appear over- luminous, i.e., located up to 0.75 mag above the sequence of single ob- 
jects. The vertical spread in the HR diagram can thus be modeled statistically to constrain 
the multiplicity frequency and the mass ratio distribution. Alternatively, individual objects 
whose spectrum contains two distinct sets of photospheric features, revealing the presence 
of a companion of much different effective temperature, are known as "spectral" binaries. 
These are easiest to detect in cases where the M* — T e s relationship is steep but the M* — 
is shallow, most notably among VLM stars and BDs. 

11 Acronyms and Definition 

• a: semi-major axis (expressed in AU unless otherwise stated) 

• BD: brown dwarf 

• CF: companion frequency (average number of companion per target in a population) 

• IMF: initial mass function 

• MF: multiplicity frequency (fraction of multiple systems in a population) 

• MS: Main Sequence 

• ONC: Orion Nebula Cluster 

• P: orbital period (expressed in days unless otherwise stated) 

• PMS: Pre-Main Sequence 

• SB: spectroscopic binary (broadly speaking, logP < 3) 

• TTS: TTauri stars 

• VB: visual binary (broadly speaking, p > 0.05") 

• VLM: very low-mass 

• a: index of the power-law fit to the orbital period distribution 

• 7: index of the power-law fit to the mass ratio distribution 

• p: angular separation of a visual binary (expressed in arcsec unless otherwise stated) 
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